Toward one-band superconductivity in MgB2 
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The two-gap model for superconductivity in MgB2 predicts that interband impurity scattering 
should be pair breaking, reducing the critical temperature. This is perhaps the only prediction of 
the model that has not been confirmed experimentally. It was previously shown theoretically that 
common substitutional impurities lead to negligible interband scattering — if the lattice is assumed 
not to distort. Here we report theoretical results showing that certain impurities can indeed produce 
lattice distortions sufficiently large to create measurable interband scattering. On this basis, we 
predict that isoelectronic codoping with Al and Na will provide a decisive test of the two-gap 
model. 

PACS numbers; 74.62.Dh,74.70.Ad,61.72.Bb,61.72. Ji 
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It is now widely accepted that MgB2 is a two-gap su- 
perconductor: its Fermi surface consists of two distinct 
sheets characterized by strong and weak electron-phonon 
couphng, respectively (see Ref.Qfor a review). This view 
is supported by numerous experiments probing either the 
larger or smaller gap, or both simultaneously. Experi- 
mental observation of the merging of the two gaps would 
constitute even stronger evidence. Such a merging is ex- 
pected, for example, from interband scattering by impu- 
rities. Within the theory of multiband superconductiv- 
ity, interband scattering mixes the "weak" and "strong" 
Cooperpairs, averaging the order parameter and reduc- 
ing Tc [31 • At small defect concentrations the suppression 
of Tc should be linear, with the larger gap decreasing and 
the smaller gap increasing. The effect should be pro- 
nounced in samples with high defect concentrations, but 
despite the relatively low quality of many samples such 
an effect has not been observed. Indeed, some samples 
with high resistivity have nearly the same critical tem- 
peratures as clean single crystals. On the other hand, 
samples doped with carbon show both gaps decreasing 
but not merging, despite a considerable reduction in Tc- 
Other types of intentional defects — whether from dop- 
ing or irradiation — have also failed to produce a merged 
gap. Even efforts to introduce defects into MgB2 for 
the explicit purpose of inducing interband scattering and 
mergingthe two gaps have failed to observe this effect 

ill in. 

This apparent lack of evidence for a central prediction 
of the two-gap model is disturbing. In Ref. [3 it was 
shown that interband scattering from substitutional im- 
purities is inherently weak, if the lattice is assumed not 
to distort. This partially explains the null results of cur- 
rent experimental efforts to induce pair breaking, but it 
does not address the possibility of pair breaking from 
impurities specifically chosen to maximize the interband 
scattering due to large lattice distortions. In this Letter 
we show that this strategy is likely to succeed and, by 
identifying a simple relationship between impurity atoms 
and the resulting lattice distortions in MgB2, suggest an 



impurity-doping protocol that will produce measurable 
pair breaking — and thereby provide the final "smoking- 
gun" evidence for the two-band model. 

Hampering such investigations is the currently limited 
insight into which defects are most effective in creating 
interband scattering. Since the states at the Fermi level 
of MgB2 are formed by the boron orbitals, one might ex- 
pect impurities (such as carbon) in the B planes to pro- 
duce large interband scattering. This is not borne out 
experimentally: substitutional C impurities have only a 
weak effect on the interband scattering jBi]. This finding 
had been anticipated theoretically as a consequence of 
the special symmetry properties of the electronic states 
within the cr band near the F point The crucial 

point is that although impurities in the B plane do have 
a strong effect on the electronic structure, they do not 
change the local point symmetry and therefore do not 
lead to significant u — tt scattering. 

The situation can be quite different for substitutions 
in the Mg plane, which may create out-of-plane distor- 
tions of B atoms in neighboring planes. Such relaxations 
do change the local point symmetry of nearby B atoms, 
mixing the in-plane Px,y and out-of-plane Pz orbitals, and 
for sufficiently large disturbances can lead to significant 
a — n scattering. Here we demonstrate by first-principles 
calculations that this is indeed the case for certain sub- 
stitutional impurities but — surprisingly — not for Mg va- 
cancies. We predict that the interband scattering effects 
will be most pronounced for isoelectronic co-doping with 
Na and Al, and that the effect on the superconducting 
properties should be detectable for impurity concentra- 
tions above 2%. 

We used density-functional theory (DFT) to study the 
lattice distortion created by Mg-plane substitutional im- 
purities from Groups I, II, and HI, by a Mg vacancy, and 
by B-plane C substitution. To model the distortion in- 
duced by single defects, we used 2x2x2 supercells of bulk 
MgB2 (and 3x3x3 supercells for convergence checks). 
Total energies and forces were calculated using projector- 
augmented- wave potentials and the generalized-gradient 
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TABLE I: First-principles displacements of nearest-neighbor 
B atoms, in angstrom, induced by various substitutional im- 
purities. Positive values indicate displacements away from the 
impurity. 



Impurity 


Bite 


In-plane, or 


(Jut-oi-plane 


Be 


Mg 


—0.014 


—0.039 


Al 


Mg 


-0.012 


-0.033 


Sc 


Mg 


+0.012 


-0.008 


C 


B 


+0.044 





Li 


Mg 


-0.006 


+0.003 


Vacancy 


Mg 


-0.005 


+0.008 


Ca 


Mg 


+0.015 


+0.028 


Na 


Mg 


+0.008 


+0.040 


K 


Mg 


+0.024 


+0.092 



approximation 0| . All atomic positions were relaxed 
within the constraint of fixed (theoretical) bulk lattice 
parameters. The resulting displacements of the nearest- 
neighbor B atoms are given in Table I. Given the strong 
intraplanar covalent bonding, it is not surprising that 
these displacements are dominated by the out-of-plane 
component (5z, which ranges from -0.04 A (for Be) to 
+0.09 A (for K); for C it is zero by symmetry. 

The out-of-plane displacements for seven different im- 
purities on the Mg site are plotted in Fig. 1. In a previous 
related study, changes in interlayer spacing induced by 
the complete substitution of Al for M g in a single plane 
were ascribed to electrostatic effects ll'. We find no 
such correlation between 5z and the formal valence of the 
impurity — for example, Be and Ca give displacements of 
opposite sign despite having identical valence — and thus 
infer that electrostatic effects are not important. On the 
other hand. Fig. 1 shows that there is an excellent cor- 
relation between 5z and the ionic radius of the impurity 
atom. Hence, we conclude that the out-of-plane displace- 
ment of B atoms by impurity atoms substituting for Mg 
is mostly a size effect. 

In light of this, one might anticipate a Mg vacancy 
to produce a large inward displacement. Our DFT re- 
sults reveal very different behavior: the vacancy creates 
a negligibly small displacement, 5z < 0.01 A. This re- 
sult is consistent with the experimental fact that pair 
breaking is not observed in low quality samples, which 
presumably contain many vacancies. However, it is very 
different from the trend shown in Fig. 1. Indeed, if one 
naively considers the vacancy as an impurity of zero size, 
the predicted displacement is almost twice that found for 
Be and Al, in sharp distinction to our DFT result. 

A simple model explains the surprisingly small dis- 
placement created by the Mg vacancy. We consider a 
MgB2 crystal containing a single substitutional defect D 
(either an impurity atom or a vacancy) in the Mg plane. 
Such a defect has 12 nearest-neighbor B atoms, consisting 
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FIG. 1: First-principles theoretical displacements of nearest- 
neighbor B atoms around a substitutional impurity in the Mg 
plane. Ionic radii are from Ref. 

of two hexagonal rings. We consider the out-of-plane dis- 
placement of the B atoms in these rings to arise from two 
opposing effects. The first of these represents the change 
in covalent bonding between B planes. In the absence of 
any defects, the MgB2 interlayer spacing cug is primar- 
ily determined by assisted hopping between B pz orbitals 
through Mg orbitals; for each B atom there are three 
such hopping paths through nearest-neighbor Mg atoms. 
With the defect present, one of these three paths now 
passes through the defect site. This new hopping path 
results in an out-of-plane force on the B atom. We as- 
sume the magnitude of this effect to be one-third of that 
found for a fully substituted Mg plane, which we approx- 
imate using the energy vs. layer spacing, E]j(c), for fully 
substituted -DB2. Thus, we consider the change in spac- 
ing between the two displaced hexagonal rings, 26z, to 
contribute an energy per B given by ^EoicMg + 'iSz). We 
have calculated the binding-energy curves Eoic) within 
DFT for AIB2, NaB2, and "vacancy-substituted" B2; the 
results are shown in Fig. 2. All can be accurately rep- 
resented by a Morse potential, which is the form we will 
use in the discussion below. 

The second effect is the restoring force experienced by 
the displaced B atoms, due to the strong covalent bond- 
ing within the B planes. For the small displacements we 
are considering it is reasonable to take this effect to be 
harmonic in Sz, again weighted by 1/3. Thus we take the 
total energy change per B to be 



1 



E{Sz) = -K{Sz 



1 



-(cMg+2(5z — cr>)/ujl2 



(1) 



where w is the width of the Morse potential, kw'^ its 
depth, and cd its equilibrium interlayer spacing. For 
small CMg — CD, it is easy to show that this energy is 
minimized for 26z = (cd — CMg)/(l + K/2k). In other 
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words, for impurities whose size is comparable to Mg, the 
displacement is linear in the size mismatch and reduced 
by the factor 1 + K/2k (which is typically in the range 
5-10); this is consistent with the DFT results shown in 
Fig. 1. Qualitatively, when cd is close to CMg, as it is 
for AIB2 and NaB2, the equilibrium displacement rep- 
resents a balance between the harmonic restoring force 



|-ftr|(5z| and a nearly harmonic (attractive or repulsive) 
force 1^ 



\k{cMg ~ Co + 2Sz) from the interlayer bonding. 
For the vacancy there is a very large size mismatch: 
Cvac is over 40% smaller than CMg- Thus for any rea- 
sonable displacement, the hexagons experience only a 
weak attractive force from the tail of the Morse potential. 
Hence, the energy is minimized for a very small displace- 
ment, Sz « —0.007 A, in agreement with the negligible 
displacement given by DFT. Qualitatively, the result of 
a large size mismatch is to largely preempt the mecha- 
nism of interlayer binding, leading to very small displace- 
ments strongly suppressed by the penalty for perturbing 
the planarity of the B layer. 

The origin of the large mismatch between CMg and Cvac 
can best be understood by comparing the band struc- 
tures of the fully substituted Z3B2 compounds at their 
equilibrium interlayer spacings. For these pure com- 
pounds (as well as the parent material) interlayer bond- 
ing arises primarily from the interaction between B pz 
orbitals in different layers. This interaction depends on 
the assisted hopping through s and Pz orbitals located on 
the D site. However, for the fully "vacancy substituted" 
compound the occupancy of the tt bands is so much re- 
duced that their contribution to bonding becomes quite 
small. At the same time, the a bands acquire substan- 
tial z-dispersion from ppn hopping, which contributes to 
bonding. Even for high vacancy concentrations (without 
complete removal of a Mg plane) it is impossible to en- 
gage the a bands in interlayer bonding by any reasonable 
dimpling of the B planes, because the planes remain too 
far apart. As a result, hardly any distortion occurs at all: 
indeed, even for 50% vacancies within a single Mg plane, 
the interlayer spacing changes by less than 0.05 A. 

We have established that the defects leading to the 
largest B displacements are Mg-plane substitutional im- 
purities with a large size mismatch (but not Mg vacan- 
cies). We now estimate the magnitude of the interband 
scattering associated with displacements from such impu- 
rities. We assume that the only relevant scattering is that 
due to the out-of-plane distortion, Sz, and use the anal- 
ogy between the formulas [l^ for the impurity-induced 
scattering rate. 
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and the electron-phonon coupling constant, 

E,.kk' <5(ek)'5(£k')l^^i',kk'lV^^iv,k-k' 
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FIG. 2: First- principles cohesive energy of completely sub- 
stituted MgB2, versus interlayer spacing. Curves are fits to 
Morse potentials. The vertical dotted line marks the equilib- 
rium interlayer spacing for MgB2. 



Here 



•-imp 



is the impurity concentration; JV^k' is the 



matrix element of the impurity perturbation potential 
(defined as the difference between the full crystal po- 
tentials with and without an impurity); ek is the elec- 
tron energy with respect to the Fermi level; w^k is the 
phonon frequency; Mk^k' is the electron-ion matrix el- 
ement (]i\dV/dQ^\k."^, where dV/dQy is the derivative 
of the crystal potential with respect to the phonon nor- 
mal coordinate Q = ^j2mLo /h,x\ and the summations are 
over all electron states and all phonon branches [l^ . 

To proceed we make three approximations, all qualita- 
tively reasonable if not quantitatively reliable. First, we 
assume that in-plane phonons, including the well-known 
E2g modes, contribute little to mterband electron- 
phonon coupling; this follows from the same symmetry 
arguments given in Ref. 0. Second, we assume by the 
same reasoning that interband impurity scattering comes 
only from the out-of-plane relaxation of B atoms. Finally, 
we approximate Eq. [3 as 

7i„,p « 12™i„p <5(£k) j {\dVlduz\^) {5zf, (4) 

and likewise Eq. |21 as 



A 



2 ^(5(£k) l\dV/duz\^) /2muo 



(5) 



Here we have assumed that the average of the crystal 
potential with respect to the vertical displacement of the 
B atom is the same in both cases. The numerical factors 
12 and 2 are the coordination of Mg and the number of 
B atoms in the unit cell, respectively. 
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First-principles calculations give the interband part of 
the electron-phonon coupling as ^^0.2 The phonon 
frequencies for the out-of-plane modes are about 400 
cm~^. Using these values, the scattering rate is given by 
7imp ~ 50ni„^p{Sz)^ eV, where Sz is in angstrom. Hence, 
for 2% Al doping we find 7inip w 1.1 meV. For 2% Na dop- 
ing we estimate a similar scattering rate, 7imp ~ 1-3 meV. 
These scattering rates are small, but should still have a 
measurable effect on the superconducting gaps and tem- 
perature. The effect on the gaps is difficult to estimate 
without full Eliashberg calculations. The reduction of 
Tc can be easily estimated using Eq. 13 from Ref. 0, 
which gives 2.0-2.5 K (in addition to any suppression 
due to the electron doping of the a band). While this 
is a small reduction compared to the changes observed 
in heavily electron-doped samples, the underlying mech- 
anism is quite different. In particular, pair breaking from 
interband scattering is unique in that it reduces Tc and 
the gap ratio while simultaneously increasing the smaller 
gap This distinctive behavior should facilitate the 
separation of interband scattering from other sources of 
Tc reduction. 

Finally, we suggest that an especially attractive test 
of these predictions would be simultaneous codoping by 
equal parts Al and Na. This would effectively be an iso- 
electronic substitution, and any effect on the supercon- 
ducting properties could then be ascribed to impurity- 
induced interband scattering. Moreover, Na and Al in- 
duce distortions of the same magnitude but of the oppo- 
site sign, which should mitigate the effects of a reduction 
in scattering due to possible short-range ordering of the 
impurities. 

In conclusion, we have performed first-principles cal- 
culations of the lattice distortion in MgB2 induced by 
several common substitutional (for Mg) impurities, and a 
Mg vacancy. We find out-of-plane displacements as large 
as 0.04 A for common impurities such as Al and Na. The 
magnitude and sign of the displacement are mainly de- 
termined by the ionic size of the impurity. On the other 



hand, for the Mg vacancy we find an essentially negligible 
displacement of nearby B atoms. The different behav- 
ior of vacancies and impurities is explained by a simple 
physical model representing the competition between in- 
terlayer binding and intralayer planarity. We estimate 
the interband scattering rate due to the Na and Al im- 
purities to be of order 1 meV, sufficiently large to give a 
detectable change in the superconducting transition tem- 
perature. Finally, we propose that the codoped mate- 
rial Mgi_2;Naa;Al2;B2, which is isoelectronic with MgB2, 
should provide an excellent test of these predictions. 
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of Naval Research. 
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